Estatística Scan e Suavização
Estatística Scan e Suavização BYM
A estatística scan foi aplicada para cada grupos de anomalias,
considerando os municípios de residência e o período de 2011-2020. É
exibido no mapa os clusters que possuem o valor p significativo. Quando
não significativo, é exibido apenas o primeiro cluster com o valor p
mais baixo.
A suavização do risco relativo foi feita através do
modelo BYM(Besag–York–Mollié), que inclui um componente ICAR para
autocorrelação espacial + um componente de efeitos aleatórios.
Defeitos de tubo neural
Scan
##
## Areas included ...........:
## Ãguas_Mornas Angelina Antônio_Carlos Biguaçu Major_Gercino Nova_Trento
## Santo_Amaro Da Imperatriz São_João Batista São_Pedro De Alcântara
## Maximum distance .........: 0.4565586
## (areas: Nova_Trento to Santo_Amaro Da Imperatriz)
## Number of cases ..........: 24
## Expected number of cases .: 10.017
## Overall relative risk ....: 2.395926
## Statistic value ..........: 7.21789
## Monte Carlo rank .........: 530
## P-value ..................: 0.053##
## Call:
## rflexscan(x = coord$x, y = coord$y, name = scan_tubo$name_muni,
## observed = scan_tubo$n, expected = scan_tubo$valor_esperado,
## nb = nb, clustersize = 10, stattype = "ORIGINAL", scanmethod = "CIRCULAR",
## simcount = 9999, rantype = "POISSON")
##
## Clusters:
## NumArea MaxDist Case Expected RR Stats P
## 1 9 0.457 24 10.017 2.396 7.218 0.0530 .
## 2 1 0.000 3 0.472 6.353 3.026 0.9644
## 3 2 0.064 7 2.556 2.739 2.632 0.9918
## 4 1 0.000 51 37.144 1.373 2.554 0.9938
## 5 2 0.121 4 1.044 3.832 2.428 0.9980
## 6 1 0.000 6 2.158 2.780 2.310 0.9988
## 7 6 0.397 8 3.447 2.321 2.206 0.9998
## 8 1 0.000 4 1.172 3.413 2.092 1.0000
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Limit length of cluster: 10
## Number of areas: 295
## Total cases: 439
## Coordinates: Cartesian
## Model: POISSON
## Scanning method: CIRCULAR
## Statistic type: ORIGINALBYM
Tabela
Microcefalia
Scan
##
## Areas included ...........:
## Catanduvas Erval_Velho Herval_D'oeste Jaborá Joaçaba Lacerdópolis
## Luzerna Ouro Salto_Veloso Treze_TÃlias
## Maximum distance .........: 0.3909241
## (areas: Ouro to Salto_Veloso)
## Number of cases ..........: 5
## Expected number of cases .: 0.9769189
## Overall relative risk ....: 5.118132
## Statistic value ..........: 4.253782
## Monte Carlo rank .........: 5393
## P-value ..................: 0.5393##
## Call:
## rflexscan(x = coord$x, y = coord$y, name = scan_microcefalia$name_muni,
## observed = scan_microcefalia$n, expected = scan_microcefalia$valor_esperado,
## nb = nb, clustersize = 10, stattype = "ORIGINAL", scanmethod = "CIRCULAR",
## simcount = 9999, rantype = "POISSON")
##
## Clusters:
## NumArea MaxDist Case Expected RR Stats P
## 1 10 0.391 5 0.977 5.118 4.254 0.5393
## 2 1 0.000 2 0.125 16.041 3.699 0.7261
## 3 3 0.210 2 0.197 10.157 2.855 0.9584
## 4 3 0.199 2 0.201 9.951 2.818 0.9634
## 5 1 0.000 1 0.023 43.459 2.801 0.9664
## 6 10 0.386 7 3.122 2.242 1.882 0.9999
## 7 1 0.000 4 1.376 2.908 1.693 1.0000
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Limit length of cluster: 10
## Number of areas: 295
## Total cases: 74
## Coordinates: Cartesian
## Model: POISSON
## Scanning method: CIRCULAR
## Statistic type: ORIGINALBYM
Tabela
Cardiopatias congênitas
Scan
##
## Areas included ...........:
## Ãguas_Mornas Antônio_Carlos Biguaçu Florianópolis
## Governador_Celso Ramos Palhoça Paulo_Lopes Santo_Amaro Da Imperatriz
## São_José São_Pedro De Alcântara
## Maximum distance .........: 0.5353028
## (areas: Florianópolis to Paulo_Lopes)
## Number of cases ..........: 264
## Expected number of cases .: 145.5059
## Overall relative risk ....: 1.81436
## Statistic value ..........: 47.17649
## Monte Carlo rank .........: 1
## P-value ..................: 1e-04##
## Areas included ...........:
## Concórdia Ipumirim Irani Lindóia_Do Sul
## Maximum distance .........: 0.3174795
## (areas: Ipumirim to Irani)
## Number of cases ..........: 37
## Expected number of cases .: 14.07358
## Overall relative risk ....: 2.62904
## Statistic value ..........: 13.10122
## Monte Carlo rank .........: 2
## P-value ..................: 2e-04##
## Call:
## rflexscan(x = coord$x, y = coord$y, name = scan_cardio$name_muni,
## observed = scan_cardio$n, expected = scan_cardio$valor_esperado,
## nb = nb, clustersize = 10, stattype = "ORIGINAL", scanmethod = "CIRCULAR",
## simcount = 9999, rantype = "POISSON")
##
## Clusters:
## NumArea MaxDist Case Expected RR Stats P
## 1 10 0.535 264 145.506 1.814 47.176 0.0001 ***
## 2 4 0.317 37 14.074 2.629 13.101 0.0002 ***
## 3 7 0.293 77 50.482 1.525 6.356 0.1267
## 4 2 0.146 5 0.596 8.395 6.243 0.1392
## 5 1 0.000 3 0.316 9.505 4.075 0.7111
## 6 1 0.000 111 86.473 1.284 3.514 0.8814
## 7 3 0.219 6 1.613 3.720 3.505 0.8837
## 8 1 0.000 4 0.831 4.812 3.121 0.9569
## 9 3 0.242 19 10.471 1.815 2.828 0.9855
## 10 9 0.527 14 8.421 1.662 1.553 1.0000
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Limit length of cluster: 10
## Number of areas: 295
## Total cases: 1022
## Coordinates: Cartesian
## Model: POISSON
## Scanning method: CIRCULAR
## Statistic type: ORIGINALBYM
Tabela
Fendas orais
Scan
##
## Areas included ...........:
## Irati Jardinópolis Serra_Alta Sul_Brasil
## Maximum distance .........: 0.1513368
## (areas: Jardinópolis to Serra_Alta)
## Number of cases ..........: 6
## Expected number of cases .: 0.9190114
## Overall relative risk ....: 6.528755
## Statistic value ..........: 6.193108
## Monte Carlo rank .........: 1410
## P-value ..................: 0.141##
## Call:
## rflexscan(x = coord$x, y = coord$y, name = scan_fendas$name_muni,
## observed = scan_fendas$n, expected = scan_fendas$valor_esperado,
## nb = nb, clustersize = 10, stattype = "ORIGINAL", scanmethod = "CIRCULAR",
## simcount = 9999, rantype = "POISSON")
##
## Clusters:
## NumArea MaxDist Case Expected RR Stats P
## 1 4 0.151 6 0.919 6.529 6.193 0.1410
## 2 6 0.323 18 7.864 2.289 4.837 0.4346
## 3 5 0.276 32 19.128 1.673 3.705 0.8233
## 4 6 0.323 11 4.744 2.319 3.020 0.9677
## 5 1 0.000 4 0.862 4.641 3.008 0.9691
## 6 3 0.145 7 2.423 2.889 2.864 0.9831
## 7 2 0.064 10 4.488 2.228 2.519 0.9970
## 8 1 0.000 2 0.304 6.575 2.073 0.9999
## 9 3 0.126 6 2.334 2.571 2.008 0.9999
## 10 7 0.640 35 24.644 1.420 1.995 0.9999
## 11 9 0.605 27 18.262 1.478 1.870 1.0000
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Limit length of cluster: 10
## Number of areas: 295
## Total cases: 771
## Coordinates: Cartesian
## Model: POISSON
## Scanning method: CIRCULAR
## Statistic type: ORIGINALBYM
Tabela
Anomalias de órgãos genitais - Hipospádias
Scan
##
## Areas included ...........:
## Ãguas_Mornas Antônio_Carlos Biguaçu Florianópolis
## Governador_Celso Ramos Palhoça Paulo_Lopes Santo_Amaro Da Imperatriz
## São_José São_Pedro De Alcântara
## Maximum distance .........: 0.5353028
## (areas: Florianópolis to Paulo_Lopes)
## Number of cases ..........: 126
## Expected number of cases .: 70.7597
## Overall relative risk ....: 1.780675
## Statistic value ..........: 21.20578
## Monte Carlo rank .........: 1
## P-value ..................: 1e-04##
## Areas included ...........:
## Alto_Bela Vista Arabutã Concórdia Ipumirim Irani Jaborá Lindóia_Do Sul
## Peritiba Presidente_Castello Branco
## Maximum distance .........: 0.4822754
## (areas: Ipumirim to Peritiba)
## Number of cases ..........: 23
## Expected number of cases .: 7.691878
## Overall relative risk ....: 2.990167
## Statistic value ..........: 10.12645
## Monte Carlo rank .........: 34
## P-value ..................: 0.0034##
## Areas included ...........:
## Chapecó Cordilheira_Alta
## Maximum distance .........: 0.08849563
## (areas: Chapecó to Cordilheira_Alta)
## Number of cases ..........: 38
## Expected number of cases .: 17.97365
## Overall relative risk ....: 2.114206
## Statistic value ..........: 8.848043
## Monte Carlo rank .........: 107
## P-value ..................: 0.0107##
## Call:
## rflexscan(x = coord$x, y = coord$y, name = scan_hipospadias$name_muni,
## observed = scan_hipospadias$n, expected = scan_hipospadias$valor_esperado,
## nb = nb, clustersize = 10, stattype = "ORIGINAL", scanmethod = "CIRCULAR",
## simcount = 9999, rantype = "POISSON")
##
## Clusters:
## NumArea MaxDist Case Expected RR Stats P
## 1 10 0.535 126 70.760 1.781 21.206 0.0001 ***
## 2 9 0.482 23 7.692 2.990 10.126 0.0034 **
## 3 2 0.088 38 17.974 2.114 8.848 0.0107 *
## 4 9 0.400 22 9.898 2.223 5.622 0.2333
## 5 1 0.000 2 0.131 15.280 3.588 0.8510
## 6 1 0.000 3 0.404 7.422 3.424 0.8924
## 7 1 0.000 3 0.564 5.319 2.584 0.9941
## 8 7 0.404 26 16.493 1.576 2.422 0.9975
## 9 1 0.000 4 1.073 3.727 2.344 0.9983
## 10 7 0.364 14 7.589 1.845 2.204 0.9994
## 11 7 0.387 9 4.178 2.154 2.108 0.9996
## 12 1 0.000 27 17.891 1.509 2.090 0.9997
## 13 1 0.000 4 1.251 3.199 1.909 1.0000
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Limit length of cluster: 10
## Number of areas: 295
## Total cases: 497
## Coordinates: Cartesian
## Model: POISSON
## Scanning method: CIRCULAR
## Statistic type: ORIGINALBYM
Tabela
Anomalias de órgãos genitais - Sexo indefinido
Scan
##
## Areas included ...........:
## Ermo Jacinto_Machado Passo_De Torres Praia_Grande Santa_Rosa Do Sul
## São_João Do Sul Turvo
## Maximum distance .........: 0.4390975
## (areas: Passo_De Torres to Turvo)
## Number of cases ..........: 5
## Expected number of cases .: 0.7945374
## Overall relative risk ....: 6.29297
## Statistic value ..........: 5.074506
## Monte Carlo rank .........: 3162
## P-value ..................: 0.3162##
## Call:
## rflexscan(x = coord$x, y = coord$y, name = scan_SI$name_muni,
## observed = scan_SI$n, expected = scan_SI$valor_esperado,
## nb = nb, clustersize = 10, stattype = "ORIGINAL", scanmethod = "CIRCULAR",
## simcount = 9999, rantype = "POISSON")
##
## Clusters:
## NumArea MaxDist Case Expected RR Stats P
## 1 7 0.439 5 0.795 6.293 5.075 0.3162
## 2 6 0.328 14 5.847 2.394 4.402 0.5077
## 3 1 0.000 1 0.021 46.888 2.873 0.9591
## 4 1 0.000 1 0.025 40.534 2.731 0.9760
## 5 1 0.000 3 0.668 4.488 2.198 0.9993
## 6 1 0.000 1 0.050 19.804 2.041 0.9998
## 7 1 0.000 1 0.060 16.586 1.873 1.0000
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Limit length of cluster: 10
## Number of areas: 295
## Total cases: 109
## Coordinates: Cartesian
## Model: POISSON
## Scanning method: CIRCULAR
## Statistic type: ORIGINALBYM
Tabela
Defeitos de membros
Scan
##
## Areas included ...........:
## Chapecó Cordilheira_Alta Guatambú
## Maximum distance .........: 0.1192868
## (areas: Cordilheira_Alta to Guatambú)
## Number of cases ..........: 127
## Expected number of cases .: 77.14058
## Overall relative risk ....: 1.646345
## Statistic value ..........: 14.08385
## Monte Carlo rank .........: 1
## P-value ..................: 1e-04##
## Areas included ...........:
## Caçador
## Maximum distance .........: 0
## (areas: Caçador to Caçador)
## Number of cases ..........: 53
## Expected number of cases .: 25.79558
## Overall relative risk ....: 2.054616
## Statistic value ..........: 11.1414
## Monte Carlo rank .........: 11
## P-value ..................: 0.0011##
## Call:
## rflexscan(x = coord$x, y = coord$y, name = scan_membros$name_muni,
## observed = scan_membros$n, expected = scan_membros$valor_esperado,
## nb = nb, clustersize = 10, stattype = "ORIGINAL", scanmethod = "CIRCULAR",
## simcount = 9999, rantype = "POISSON")
##
## Clusters:
## NumArea MaxDist Case Expected RR Stats P
## 1 3 0.119 127 77.141 1.646 14.084 0.0001 ***
## 2 1 0.000 53 25.796 2.055 11.141 0.0011 **
## 3 7 0.404 100 68.960 1.450 6.365 0.1264
## 4 6 0.290 32 17.714 1.806 4.688 0.5173
## 5 1 0.000 6 1.253 4.789 4.657 0.5294
## 6 1 0.000 15 6.266 2.394 4.378 0.6240
## 7 9 0.494 96 71.321 1.346 4.001 0.7536
## 8 10 0.535 340 295.852 1.149 3.693 0.8510
## 9 6 0.282 23 12.811 1.795 3.296 0.9404
## 10 7 0.272 38 24.734 1.536 3.094 0.9685
## 11 2 0.093 11 4.765 2.309 2.977 0.9786
## 12 1 0.000 3 0.664 4.520 2.191 0.9999
## 13 1 0.000 3 0.745 4.027 1.925 1.0000
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Limit length of cluster: 10
## Number of areas: 295
## Total cases: 2078
## Coordinates: Cartesian
## Model: POISSON
## Scanning method: CIRCULAR
## Statistic type: ORIGINALBYM
Tabela
Defeitos da parede abdominal
Scan
##
## Areas included ...........:
## Dona_Emma Doutor_Pedrinho José_Boiteux Vitor_Meireles Witmarsum
## Maximum distance .........: 0.3915501
## (areas: Doutor_Pedrinho to Witmarsum)
## Number of cases ..........: 6
## Expected number of cases .: 1.107336
## Overall relative risk ....: 5.418411
## Statistic value ..........: 5.282066
## Monte Carlo rank .........: 2821
## P-value ..................: 0.2821##
## Call:
## rflexscan(x = coord$x, y = coord$y, name = scan_abdominal$name_muni,
## observed = scan_abdominal$n, expected = scan_abdominal$valor_esperado,
## nb = nb, clustersize = 10, stattype = "ORIGINAL", scanmethod = "CIRCULAR",
## simcount = 9999, rantype = "POISSON")
##
## Clusters:
## NumArea MaxDist Case Expected RR Stats P
## 1 5 0.392 6 1.107 5.418 5.282 0.2821
## 2 1 0.000 2 0.141 14.141 3.445 0.8808
## 3 5 0.305 8 2.729 2.932 3.375 0.8959
## 4 7 0.382 13 5.847 2.224 3.313 0.9075
## 5 1 0.000 3 0.471 6.366 3.034 0.9583
## 6 10 0.401 12 5.527 2.171 2.894 0.9739
## 7 3 0.134 3 0.578 5.192 2.528 0.9960
## 8 1 0.000 3 0.791 3.791 1.796 1.0000
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Limit length of cluster: 10
## Number of areas: 295
## Total cases: 336
## Coordinates: Cartesian
## Model: POISSON
## Scanning method: CIRCULAR
## Statistic type: ORIGINALBYM
Tabela
Síndrome de Down
Scan
##
## Areas included ...........:
## Catanduvas Concórdia Irani Jaborá Joaçaba Lacerdópolis Ouro
## Presidente_Castello Branco
## Maximum distance .........: 0.3740191
## (areas: Catanduvas to Concórdia)
## Number of cases ..........: 27
## Expected number of cases .: 10.6338
## Overall relative risk ....: 2.539074
## Statistic value ..........: 9.045609
## Monte Carlo rank .........: 99
## P-value ..................: 0.0099##
## Call:
## rflexscan(x = coord$x, y = coord$y, name = scan_down$name_muni,
## observed = scan_down$n, expected = scan_down$valor_esperado,
## nb = nb, clustersize = 10, stattype = "ORIGINAL", scanmethod = "CIRCULAR",
## simcount = 9999, rantype = "POISSON")
##
## Clusters:
## NumArea MaxDist Case Expected RR Stats P
## 1 8 0.374 27 10.634 2.539 9.046 0.0099 **
## 2 10 0.403 39 23.329 1.672 4.607 0.4972
## 3 2 0.004 5 0.948 5.273 4.276 0.6134
## 4 3 0.196 55 37.021 1.486 4.115 0.6727
## 5 1 0.000 15 6.524 2.299 4.079 0.6872
## 6 1 0.000 3 0.329 9.131 3.970 0.7274
## 7 2 0.077 4 0.647 6.185 3.945 0.7364
## 8 7 0.498 27 15.516 1.740 3.598 0.8456
## 9 1 0.000 2 0.166 12.048 3.147 0.9459
## 10 1 0.000 2 0.185 10.809 2.949 0.9734
## 11 1 0.000 27 18.376 1.469 1.836 1.0000
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Limit length of cluster: 10
## Number of areas: 295
## Total cases: 545
## Coordinates: Cartesian
## Model: POISSON
## Scanning method: CIRCULAR
## Statistic type: ORIGINAL